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ABSTRACT 

We present the AGILE-GRID monitoring of Cygnus X-3, during the period between November 2007 and July 2009. We report here the whole 
AGILE-GRID monitoring of Cygnus X-3 in the AGILE "pointing" mode data-taking, in order to confirm that the y-ray activity occurs in coinci- 
dence with the same repetitive pattern of multiwavelength emission and to analyze in depth the overall y-ray spectrum by assuming both leptonic 
and hadronic scenarios. Seven intense y-ray events were detected in this period, with a typical event lasting 1 or 2 days. Such a duration is longer 
than the likely cooling times of the y-ray emitting particles, implying we are seeing continuous acceleration rather than the result of an impulsive 
event such as the ejection of a single plasmoid which then cools as it propagates outwards. Cross-correlating the AGILE-GRID light curve with 
X-ray and radio monitoring data, we find that the main events of y-ray activity have been detected while the system was in soft spectral X-ray states 
(RXTE/ ASM count rate > 3 counts s _1 ), coincident with local and often sharp minima of the hard X-ray flux (SwiftfB AT count rate < 0.02 counts 
crrT 2 s~'), a few days before intense radio outbursts. This temporal repetitive coincidence between y-ray transient emission and spectral state 
changes of the source turns out to be the spectral signature of the y-ray activity from this microquasar. The y-ray events thus may reflect a sharp 
transition in the structure of the accretion disk and its corona, which leads to a rebirth of the microquasar jet and subsequent enhanced radio 
activity. The y-ray differential spectrum of Cygnus X-3 (100 MeV - 3 GeV), obtained by averaging the data collected by the AGILE-GRID during 
the y-ray events, is consistent with a power law of photon index a = 2.0 ± 0.2. Finally, we examine leptonic and hadronic emission models for the 
y-ray events and find that both scenarios may work. In the leptonic model - based on IC scatterings of mildly relativistic electrons on soft photons 
from the Wolf-Rayet companion star and from the accretion disk - the emitting particles may also contribute to the overall hard X-ray spectrum, 
possibly explaining the hard non- thermal power-law tail sometimes seen during special soft X-ray states in Cygnus X-3. 

Key words, stars: individual: Cygnus X-3 - gamma rays: general - X-rays: binaries - radio continuum: general - radiation mechanisms: non- 
thermal - stars: winds, outflows 



1. Introduction 

Cygnus X-3 is the brightest radio source among microquasars 
and it was discovered, as an X-ray source, in 1966 (Giacconi et 
al., 1967). It is a high mass X-ray binary: the companion star is a 
Wolf-Rayet (WR) star (van Kerkwijk et al., 1992) with a strong 
stellar helium wind (Szostek & Zdziarski , 2008). The system 
is located at a distance of about 7-10 kpc (Bonnet-Bidaud & 
Chardin, 1988; Ling et al., 2009). The orbital period is 4.8 hours, 
inferred from infrared (Becklin et al., 1973), X-ray (Parsignault 
et al., 1972) and y-ray (Abdo et al., 2009) observations. Due 
to the very tight orbit (orbital distance d ~ 3 x 10 11 cm), 
the compact object is totally enshrouded in the wind of the 
companion star 1 . The nature of the compact object is still 
uncertain 2 (Vilhu et al., 2009), even if a black hole scenario is 
favored (Szostek & Zdziarski , 2008; Szostek et al., 2008). In 



1 The observational evidence of this strong wind can be found in 
the prominent attenuation of the Cygnus X-3 power density spectrum 
(PDS) for frequencies above 0.1 Hz (Axelsson et al., 2009; Koljonen et 
al., 2011). 

2 Published results suggest either a neutron star of 1.4 M Q (Stark & 
Saia, 2003) or a black hole with a mass < 10 M Q (Hanson et al., 2000; 
Shrader et al., 2010). 



the radio band the system shows strong flares ("major radio 
flares") reaching up to few tens of Jy. Radio observations at 
milliarcsec scales confirm emissions (at cm wavelengths) from 
both a core and a one-sided relativistic jet (v ~ 0.81c), with an 
inclination to the line-of-sight of < 14° (Mioduszewski et al., 
2001). The radiation from the jet dominates the radio emission 
from the core during (and soon after) the major flares (Tudose 
etal., 2010). 

Cygnus X-3 exhibits a clear and repetitive pattern of 
(anti)correlations between radio and X-ray emission, and 
an overall anticorrelation between soft and hard X-ray fluxes 
(McCollough et al., 1999; Szostek et al., 2008). The most 
important pattern of correlations found by Szostek et al. (2008) 
is related to the connection between radio (8.3 GHz band, GBI) 
and soft X-ray emissions (3-5 keV band, RXTE/ASM). When 
the soft X-ray flux is above the transition level (3 counts/s), the 
source can be found in different states, depending on the level 
of the radio flux density. In particular, the quenched state is 
characterized by a radio flux density < 30 mJy and it is followed 
by a major-flaring state with values of the radio flux density 
> 1 Jy. It is very important to emphasize that all major radio 
flares have been observed after a quenched state, and in almost 
all cases the quenched state is followed by a major flare. After a 
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major flare a "hysteresis" in the radio/soft-X-ray plane is found: 
the decline of the radio flux density never occurs through a 
quenched state. 

Firm detections of high-energy y-rays (HE y-rays: >100 
MeV) from Cygnus X-3 3 have been published at the end of 
2009: the AGILE Team found evidence of strong y-ray tran- 
sient emission above 100 MeV in coincidence with special X- 
ray/radio spectral states (Tavani et al., 2009a), and the Fermi- 
LAT collaboration announced the detection of the y-ray orbital 
modulation (Abdo et al., 2009). The peak y-ray isotropic lumi- 
nosity detected above 100 MeV is L y ~ 10 36 erg s (for a dis- 
tance of 7-10 kpc). The y-ray emission is most likely associ- 
ated with relativistic jet (Tavani et al., 2009a; Abdo et al., 2009; 
Dubus et al., 2010; Cerutti et al., 2011; Zdziarski et al., 2012), 
but the radiative process (leptonic or hadronic) is not certain. 

A possible leptonic scenario for y-ray emission in 
Cygnus X-3 was proposed by Dubus et al., 2010: stellar UV pho- 
tons are Compton upscattered to HE by relativistic electrons ac- 
celerated in the jet. The particle acceleration could take place 
in a shock where the jet interacts with the dense stellar wind of 
the WR star. The jet - with moderate bulk relativistic speed - 
is oriented not too far from the line-of-sight. According to this 
model, Cygnus X-3 is a microblazar, with a jet pointing towards 
the Earth. 

The y-ray modulation - coherent with the orbital period - sug- 
gests that the emitting region is located between ~10 10 cm and 
~3 x 10 12 cm (10c/) to the compact object (Dubus et al., 2010; 
Cerutti et al., 201 1). The lack of modulation at radio wavelengths 
and the delay (~5 days, Abdo et al., 2009) between the onset of 
y-ray activity and the radio flare suggest different emission re- 
gions linked by the collimated jet. The y-ray emission - related 
to inverse Compton (IC) scatterings - most likely occurs close 
to the compact object, while the radio emission - assumed to be 
synchrotron in origin - occurs farther out in the jet, at an angular 
distance from the core of few tens of milli-arcsec (e.g., Tudose 
et al., 2007, 2010), corresponding to ~ 10 15 -10 16 cm. The y-ray 
modulation is due to the anisotropic efficiency of the IC scat- 
tering (Aharonian & Atoyan, 1981). Thus, the y-ray maximum 
occurs at the superior conjuction (compact object behind the WR 
star), when relativistic electrons of the jet, moving towards the 
Earth, have head-on collision with UV stellar photons. This or- 
bital phase corresponds to the minimum of the X-ray modula- 
tion, due to the maximum absorption/scattering by the compan- 
ion's wind (Abdo et al., 2009; Dubus et al., 2010; Zdziarski et 
al., 2012). 

A hadronic scenario accounting for y-ray emission in the 
context of microquasars was discussed by e.g. Romero et al. 
(2003, 2005). The model is based on the interaction of a mildly 
relativistic jet with the dense wind of the companion star, y-ray 
emission is due to the decay of neutral pions (n ) produced by 
pp collisions. 



3 y-ray detections of Cygnus X-3 were reported in the 1970s and in 
the 1980s at TeV (Vladimirsky et al., 1973; Danaher et al., 1981; Lamb 
et al., 1982) and PeV energies (Samorski & Stamm, 1983; Bhat et al., 
1986). However, subsequent observations by more sensitive ground- 
based telescopes did not confirm TeV and PeV emission from the 
source (O'Flaherty et al., 1992). Furthermore, the COS-B satellite could 
not find any clear emission from Cygnus X-3 at MeV-GeV energies 
(Hermsen et al., 1987), and both CGtfO/EGRET observations of the 
Cygnus region (1991-1994) and the first-year analysis of AGILE obser- 
vations could not demonstrate a solid association with the microquasar, 
although confirming a y-ray detection above 100 MeV in a region in- 
cluding Cygnus X-3 (Mori et al., 1997; Pittori et al., 2009). 



Furthermore, TeV emission from relativistic jet in micro- 
quasars has been predicted by several models (e.g., see Atoyan & 
Aharonian, 1999). A search for very-high-energy (VHE) y-rays 
from microquasar GRS 1915+105 with the HESS Cherenkov 
telescope has been carried out, but no significant detection was 
found in the direction of the source (H.E.S.S. Collaboration et 
al., 2009). On the other hand, hints of VHE y-rays were found in 
Cygnus X-l (Albert et al., 2007). MAGIC observed Cygnus X-3 
several times between March 2006 and August 2009, during 
both hard and soft states , but no evidence of clear VHE y-ray 
emission from the microquasar was found: an overall 2cr upper 
limit to the integral flux has been set to 2.2 x 10~ 12 photons cm -2 
s _1 for energies above 250 GeV (Aleksic et al., 2010). 

Here we present a comprehensive and homogeneous analysis 
of Cygnus X-3 which takes into account y-ray events found in 
the data between 2007-November-02 and 2009-July-29, during 
the AGILE "pointing" mode data-taking. We analyzed a dataset 
previously published by Tavani et al. (2009a) and Bulgarelli et al. 
(2012a). We report here the whole AGILE-GRID monitoring of 
Cygnus X-3 during the "pointing" mode, in order to confirm that 
the y-ray activity occurs in coincidence with the same repetitive 
pattern of multiwavelength emission and to analyze in depth the 
overall y-ray spectrum by assuming both leptonic and hadronic 
scenarios. 

2. Observations 

The AGILE scientific instrument (Tavani et al., 2009b) is very 
compact and is characterized by two co-aligned imaging detec- 
tors operating in the energy ranges 30 MeV-30 GeV (GRID: 
Barbiellini et al., 2002; Prest et al., 2003) and 18-60 keV (Super- 
AGILE: Feroci et al., 2007), as well as by an anticoincidence 
system (Perotti et al., 2006) and a calorimeter (Labanti et al., 
2006). AGILE'S performance is characterized by large fields of 
view (2.5 and 1 sr for the y-ray and hard X-ray bands, respec- 
tively) and optimal angular resolution (PSF = 3.0° at 100 MeV 
and PSF = 1.5° at 400 MeV, see Cattaneo et al., 201 1). 
AGILE operated in "pointing" mode with fixed attitude until 
mid-October, 2009; in November, 2009 AGILE entered in "scan- 
ning mode", characterized by a controlled rotation of the point- 
ing axis. 

During the "pointing" mode data-taking (~2.5 years) the AGILE 
satellite performed ~100 pointings with variable exposure times 
(typically 3-30 days), drifting about 1 degree per day with re- 
spect to the initial boresight direction to match solar panel il- 
lumination constraints 5 . In this configuration the AGILE-GRID 
was characterized by enhanced performances in the monitoring 
capability of a given source, especially in the energy band 100- 
400 MeV (see Bulgarelli et al., 2012a for details). Due to the 
different pointing strategies between the AGILE and Fermi satel- 
lites, the high on-source cumulative exposure (between 100 and 
400 MeV) of the AG]LE-GRID may be fundamental in the ob- 
servation of this particular source. 

In this paper we report an analysis based on the AGILE- 
GRID data collected between 2007-November-02 and 2009- 
July-29 (the same dataset reported by Chen, Piano et al., 201 1). 

4 The MAGIC telescope also pointed at Cygnus X-3 after two y-ray 
alerts from the AGILE-GRID team (the first one after the y-ray event 
of 16-17 April 2008, the second after the event of 13-14 July 2009, see 
Appendix A). In both cases they found a 2<x upper limit, for energies 
above 250 GeV, of ~10~ n photons cnr 2 s" 1 (Aleksic et al., 2010). 

5 A detailed schedule of the AGILE observations - with the pointing 
starting coordinates and observation's starting and ending times - is 
available online at http://agile.asdc. asi.it/current_pointing.html 
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During this period AGILE repeatedly pointed at the Cygnus re- 
gion for a total of ~275 days, corresponding to a net exposure 
time of ~1 1 Ms. The detailed analysis of the dataset is presented 
in Appendix A. In this paper we report 7 y-ray main events. 



3. The gamma-ray activity in a multiwavelength 
context 

3.1. General characteristics of the gamma-ray events 

Figure 1 shows the comprehensive multiwavelength light curve 
of Cygnus X-3, in order to analyze the pattern of emission. The 
y-ray activity detected by the AGILE-GRID is presented along 
with hard X-ray fluxes from Swift/BAT 6 (15-50 keV), soft X- 
ray fluxes from RXTE/ASM 1 (3-5 keV) and radio flux density 
(when available) from AMI-LA 8 (15 GHz) and RATAN-600 9 
(2.15, 4.8, 11.2 GHz) radio telescopes. 

The aim of the cross-correlation between the AGILE-GRID 
light curve and the multiwavelength emission pattern is to dis- 
cuss the y-ray trigger criteria and to compare them with the ones 
previously published by Tavani et al., 2009a; Bulgarelli et al., 
2012a; Corbel et al., 2012. 

Observing the light curve in Figure 1, as well as the detailed 
zooms in Figure 2, we notice that: 

- there is a strong anticorrelation between hard X-ray and 
y-ray emission. Every local minimum of the hard X-ray 
light curve is associated with y-ray emission detected by 
the AGILE-GRID (see also the weak y-ray event detected 
on 2008-June-21, MJD = 54638.58, in the plot of Figure 1: 
yffS =2.77, photon flux = (131 + 61) x 10~ 8 photons cm" 2 
s~'). Vice versa, every time the AGILE-GRID detects y-ray 
activity the system exhibits a deep local minimum of the 
hard X-ray light curve (Swift/BAT count rate < 0.02 counts 
cirr 2 s _1 ). 

- Every time we detect y-ray activity, Cygnus X-3 is in a soft 
spectral state (RXTE/ASM count rate > 3 counts s _1 , i.e., the 
transitional level defined by Szostek et al., 2008). 

- Every time we detect y-ray episodes (Table A.l and red 
points in the AGTLE-GRID light curve in Figure 1) the sys- 
tem is moving to a major radio flare (radio flux density > 1 
Jy) or to a quenched state preceding a major radio flare. 

In Table 1 we report a brief synthesis of the multi-frequency pat- 
tern of emission of each main y-ray event ( VrS > 3) detected 
by the AGTLE-GRID. The average value of the delay between 
the y-ray event and the subsequent radio flare is AT2 ~ 4 days, 
consistent with the value (5+7 days) found by Abdo et al. (2009). 
If we refer to the third column in Table 1, we can see that the 
transient y-ray emission occurs when the system is moving into a 
quenched state ("Pre-Quenched") or towards a radio flare ("Pre- 
Flare") that has always been observed after a quenched state: i.e., 
the y-ray emission is detected when the system is moving into 
or out of a quenched state. So that, from a purely phenomeno- 
logical point of view, the quenched state seems to be a "key" 
condition for the y-ray emission. 



6 Swift/BAT (Burst Alert Telescope) transient monitor results pro- 
vided by the Swift/BAT team. 

7 Rossi X-ray Timing Explorer (RXTE), All-Sky Monitor (ASM). 
Quick-look results provided by the RXTE/ASM team. 

8 Courtesy of the Arcminute Microkelvin Imager (AMI) Team. 

9 Courtesy of S. Trushkin and the RATAN-600 Team. 



Table 1: Main events of y-ray activity ( y/TS > 3) detected by the 
AGILE-GRID. Column one: date of the y-ray event (average in MJD). 
Column two: X-ray spectral state. Column three: radio-flux-density state 
at the time of the y-ray activity. Column four: time delay (AT2) in days 
between the y-ray event and the major radio flare. Column five: radio 
flux density of the major radio flare. 



MJD 


X-ray 
State 


Radio 
State 


AT 2 
[days] 


Radio 
flux density 


54507.19 


X-ray state 
trans, level 


(?) 


(?) 


(?) 


54573.08 


Soft 


Pre-Flare 


~ 1 


«16Jy (11.2 GHz) 


54773.17 


Soft 


Pre-Quenched 


~ 8 


KlJy(15GHz) 


54812.39 


Soft 


Pre-Flare 


~ 6 


«3 Jy (1 1.2 GHz) 


55003.37 


Soft 


Pre-Quenched 


(?) 


(?) 


55025.55 


Soft 


Pre-Flare 


~ 3 


«3 Jy (11.2 GHz) 


55034.88 


Soft 


Pre-Flare 


~ 1 


«2.3 Jy (11.2 GHz) 
»1.6Jy(15GHz) 



We stress that in general y-ray events - always during soft 
states - occur in proximity of spectral X-ray transitions. In par- 
ticular we found that intense transient y-ray emission is de- 
tected both immediately after hard-to-soft (e.g., the events of 2-3 
November 2008 [MJD = 54773.17] and 20-21 June 2009 [MJD 
= 55003.37]) and before soft-to-hard spectral X-ray transitions 
(e.g., the events of 16-17 April 2008 [MJD = 54573.08], 11-12 
December 2008 [MJD = 54812.39] and 13-14 July 2009 [MJD 
= 55025.55]). Observing Figure 1 and the zooms in Figure 2, 
we can note that strong transient y-ray emission generally oc- 
curs when the system has just entered into or is moving out 
of a prominent minimum of the Swift/BAT light curve (e.g., 
November-December 2008 [MJD * 54770-54815] and June- 
July 2009 [MJD =s 55000-55040] events.). 

This comprehensive study confirms that the y-ray emission 
conditions of Cygnus X-3, during the whole "pointing" moni- 
toring by the AGILE satellite, are completely in agreement with 
the ones found by Tavani et al., 2009a; Bulgarelli et al., 2012a; 
Corbel et al., 2012. 

3.2. Peculiarities of the gamma-ray events 

3.2.1. The gamma-ray event of 1 1-12 February 2008 
(MJD = 54507.19) 

A more detailed discussion is needed for the y-ray event of 1 1- 
12 February 2008 (MJD = 54507.19, see the upper-left panel 
in Figure 2), which is a special event among the AGTLE-GRID 
detections. 

The y-ray activity in this period occurred when Cygnus X-3 
was at the transitional level between soft and hard X-ray spec- 
tral state (RXTE/ASM count rate ~3 counts s~') and coincides 
with a a short but sharp dip of the Swift/BAT light curve: MJD = 
54508, count rate = (0.019 +0.005) counts cirr 2 s" 1 , DATA_FLAG 
= (data quality flag = good). Even if the event does not occur 
in coincidence with a bright soft state, it seems to confirm the 
simultaneous y-ray-event/hard-X-ray-minimum occurrence that 
we find in all other cases. There are no available radio data for 
this period. 
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Fig. 1: Multi-frequency light curve of Cygnus X-3 from 2007-December-12 to 2009-September-26 (MJD: 54450-55100). From top to bottom: radio flux density [RATAN-600 (2.15, 4.8, 1 1.2 GHz) 
and AMI-LA (15 GHz)], soft X-ray count rate [RXTE/ASM (3-5 keV)], hard X-ray count rate [Swift/BAT (15-50 keV)] and gamma-ray photon fluxes [AGILE-Gfi/£> (above 100 MeV)]. In the 
bottom panel gray regions represent the AGILE pointing at the Cygnus region; magenta points are the main events of y-ray activity with yTS > 3 (see Table A. 1), black points are the y-ray 
detections with 2 < y/TS < 3 and dark-gray arrows are the 2cr upper limits related to yTS < 2. The dashed lines in the panels of the RXTE/ASM and Swift/BAT count rate represent the transition 
level of 3 counts s -1 and 0.02 counts cirr 2 s _I respectively (see Section 3.1 for details). 
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Fig. 2: Multi-frequency light curves centered on the main events of y-ray activity detected by the AGILE-GRID (detailed views of the main plot in 
Figure 1). Upper-left plot: from 2008-January-30 to 2008-February-24 (MJD: 54495-54520). Upper-right plot: from 2008-April-04 to 2008-May- 
04 (MJD: 54560-54590). Lower-left plot: from 2008-October-ll to 2009-January-09 (MJD: 54750-54840). Lower-right plot: from 2009-June-08 
to 2009-August-02 (MJD: 54990-55045). 
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3.2.2. The gamma-ray event of 20-21 June 2009 
(MJD = 55003.37) 

The y-ray event on 20-21 June 2009 (MJD = 55003.37) occurred 
at the beginning of a quenched radio state (see the lower-right 
panel in Figure 2); unfortunately we have no subsequent radio 
data covering the following days (from the 3rd to the 9th day 
after the y-ray event), but we cannot exclude the presence of 
a major radio flare - soon after the quenched state - that might 
have escaped detections. In Figure 2 (Lower-right panel) we note 
that the AMI-LA light curve (15 GHz) shows a following radio 
flare ~ 1 8 days after the y-ray event (~ 1 3 days after the quenched 
state), possibly indicating a delayed radio burst (MJD = 55021 .2, 
radio flux density = 1.06 Jy). 

3.2.3. The long-term gamma-ray emission June-July 2009 
(MJD * 55000-55040) 

Finally, we note that the extended y-ray activity of June-July 
2009 (see the lower-right panel in Figure 2), occurred during an 
intense radio activity (with a high average radio flux density) 
coincident with a long-lasting soft X-ray spectral state. It is im- 
portant to remark that this period coincides with one of the two 
temporal windows of strong y-ray activity detected by Fermi- 
LAT, with a peak photon flux greater than ~ 200 x 10~ 8 photons 
ctrT 2 s~' (Abdo et al., 2009). This y-ray peak is simultaneous 
(and consistent) with the event detected by the AGILE-GRID on 
21-23 July 2009 (photon flux = (158 + 59) x 10~ 8 photons ctrT 2 
s~'). These data (June-July 2009) are discussed in more detail 
in a dedicated paper on the y-ray emission of the microquasar 
during the mid-2009/mid-2010 period (Bulgarelli et al., 2012a). 



4. Modeling the Spectral Energy Distribution 

By accounting for the X-ray, y-ray (AGILE-GRID) and TeV 
emission (MAGIC spectral upper limits), we modeled the multi- 
wavelength Spectral Energy Distribution (SED) of Cygnus X-3 
during soft spectral state, with both leptonic and hadronic sce- 
narios. We considered a X-ray spectrum measured by RXTE- 
PCA 10 and RXTE-HEXTE U (-3-150 keV) when the source 
was in a "hypersoft" state (Koljonen et al., 2010), the AGILE- 
GRID spectrum for the main y-ray events (Figure A.l), and 
the MAGIC differential flux upper limits obtained when the 
source was in the soft state. (Aleksic et al., 2010). The hyper- 
soft state of Cygnus X-3, subclass of the ultrasoft state defined 
in Hjalmarsdotter et al., 2009, is usually exhibited by the mi- 
croquasar during the quenching/pre-flaring radio activity 12 . This 
X-ray spectral state is characterized by a weak and hard power 
law tail (a- 1.7-1.9) with a non-thermal origin. 
Thus, we analyzed a pattern of multiwavelength datasets that, 
even if not simultaneous, are qualitatively consistent because 
they are all referred to the same spectral state of Cygnus X-3: 
the X-ray and AGILE-GRID datasets are related to the soft-state 
activity preceding the radio major flares, the TeV data are re- 
lated to the soft-state activity following the radio major flares 
(MAGIC has never observed Cygnus X-3 during the pre-flaring 
radio states). 



10 Proportional Counter Array (PCA) 

1 1 High Energy X-ray Timing Experiment (HEXTE) 

12 This is an average "hypersoft" spectrum related to 28 pointed 
RXTE observation between February 2000 and January 2006 (see the 
Supporting Information of Koljonen et al., 2010 for details). 



4.1. A leptonic scenario 

We have modeled the multi-frequency SED by assuming a 
simple leptonic scenario: a plasmoid of high energy elec- 
trons/positrons, injected in the jet structure, that upscatters - via 
inverse Compton interactions - soft seed photons from the WR 
star and from the accretion disk. 

Our aim is to analyze a possible link between the power law 
tail of this special soft X-ray spectral state and the y-ray emission 
detected by the AGILE-GRID. 

The physical parameters of the photon field are literature- 
based. We modeled the X-ray data with a black body (BB) spec- 
trum 13 characterized by a a; 1 .3 keV (consistent with the typ- 
ical value of the characteristic temperature of the disk during the 
hypersoft/ultrasoft state: Hjalmarsdotter et al., 2009; Koljonen 
et al., 2010) and a L&j ~ 8 x 10 37 erg s _1 . The main parameters 
we used for the WR star are: T+ = 10 5 K, oi 10 39 erg s _1 
(see Dubus et al., 2010). The WR star is assumed to emit UV 
photons isotropically. We modeled the average y-ray emission 
in the orbital phase, so we assumed the seed photons from the 
companion star colliding with the leptons mainly on the side of 
the jet. 

We have carried out two different models: in the first one 
(leptonic model "A") the plasmoid interacts with the soft pho- 
ton bath "close" to the disk (the star-plasmoid distance is R « 
d ss 3 x 10 11 cm), in the second one (leptonic model "B") the 
interaction region is "far away" from the accretion disk (the star- 
plasmoid distance is R « 10a! ~ 3 x 10 12 cm). 
For both models, the inclination of the jet to the line of sight is 
assumed to be i = 14°, and the plasmoid is assumed to be spheri- 
cal (radius r = 3 xlO 10 cm) with a bulk motion characterized by a 
Lorentz factor of T — 1.5 (v = V5c/3). The population of elec- 
trons is modeled by a broken-power-law spectral distribution, 
with spectral indices a\ = 2.2, ai = 4.0, y m ,„ = 1, y max = 10 5 
and energy break of y^ = 4 x 10 3 : 



dN 
dydV 



K, 



\y,„ 



J (1) 



The spectral indices and the energy break of the electron distri- 
bution are the best-fitting values on the AGILE-GRID spectral 
shape. The distribution of electrons/positrons is assumed to be 
isotropic in the plasmoid rest frame (the jet comoving frame). 
We adopted the Klein-Nishina formula to describe the Compton 
scattering of soft photons by a cloud of mildly relativistic leptons 
(Aharonian & Atoyan, 1981). 

In the leptonic model "A" the distance from the star to the 
plasmoid location is assumed to be R as 3 x 10 11 cm (R as d), 
i.e., the plasmoid in the jet is very close to the compact object 
and to the accretion disk. The distance between the plasmoid 
center and the compact object is H as 3 x 10 10 cm, i.e., H « r. 
The result of this modeling is in Figure 3. The electron num- 
ber density of the plasmoid is n e as 3 x 10 9 electrons ctrr 3 (the 
prefactor K e in Eq. (1) is K e — 2 x 10 5 ctrr 3 , the integrated num- 
ber of electrons is N e — 3 x 10 41 ). We took into account the yy 
absorption (for e ± pair production) between the X-ray photons 



13 Our modeling is a simplification: we are assuming that the bump in 
the X-ray emission during the "hypersoft" state can be modeled with a 
simple BB component, that is a very good approximation for our pur- 
pose. In this state the overall X-ray emission is totally dominated by a 
strong BB emission from the accretion disk. Nevertheless, a more accu- 
rate modeling should be based on a Comptonized BB spectrum by the 
corona (see Koljonen et al., 2010 for details). 



6 



G. Piano et al.: The AGILE monitoring of Cygnus X-3: transient gamma-ray emission and spectral constraints 



from disk and the IC y-rays photons. We assumed that the dis- 
tribution of the disk photons is fully isotropized by the stellar 
wind in the observer frame. This implies that the y-ray photo- 
sphere (i.e., where r yy > 1) has a radius of ~ 10 10 cm (Cerutti et 
al., 2011). With these assumptions, the lowest part of the plas- 
moid is within the y-ray photosphere 14 . The spectral component 
related to IC scatterings on the disk photons (green curve) is ac- 
tually produced in this region, very close to the disk, where the 
X-ray photon density as well as the optical depth is high. Since 
T yy > 1, this component shows a sharp cut-off at ~100 MeV (i.e., 
the threshold for e ± production, given the characteristic energies 
of the disk photons). On the other hand, the spectral component 
related to IC scatterings on the stellar wind photons (red curve) 
does not show any cut-off, because it is mainly produced in the 
farthest part of the plasmoid (outside the y-ray photosphere, for 
distances greater than ~ 10 10 cm with respect to the disk), where 
the yy absorption by the X-ray disk photons is negligible. Thus, 
we deduce that in our geometry the plasmoid volume outside the 
y-ray photosphere emits the bulk of the y-ray emission above 
100 MeV via IC processes on stellar photons (see Figure 3). 
In model "A", assuming a lepton injection rate of N e = n e nr 2 v as 
2 x 10 41 leptons s , the jet kinetic luminosity for the leptons 
(L kin> e = N e T m e c 2 ) would be L^ n e as 2 x 10 35 erg s 4 . 

In the leptonic model "B " the distance from the star to the 
plasmoid is assumed to be R w 3 x 10 12 cm (R as lOd), i.e., the 
plasmoid in the jet is far away from the compact object and the 
accretion disk. The distance between the plasmoid center and the 
compact object is H at 3 x 10 12 cm, i.e., R ~ H. We assumed 
that the disk photons are entering in the plasmoid mainly from 
behind. The results of this modeling are shown in Figure 4. The 
electron density of the plasmoid is n e as 1.5 x 10 11 electrons 
cm -3 (the prefactor K e in Eq. (1) is K e - 8 x 10 6 cm 4 , the in- 
tegrated number of electrons is N e - 1.5 x 10 43 ). In this model 
the spectral component related to the IC scatterings on disk pho- 
tons (green curve) is negligible compared to the IC component 
on soft photons from the star (red curve). We remark that the "IC 
disk" component does not show any cut-off related to yy absorp- 
tion by X-ray photons, because the IC y-rays are produced well 
outside the y-ray photosphere (at distances » 10 10 cm). 
In model "B", assuming a lepton injection rate of N e = 
n e nr 2 v as 10 43 leptons s _1 , the jet kinetic luminosity for the lep- 
tons (L kin , e=N e Y m e c 2 ) would be e as 10 37 erg s _I . 

In these models, the expected VHE y-ray emission would be 
very faint. These expectations are consistent with the MAGIC 
upper limits, and might explain the lack of TeV bright detections 
during soft states. 

4.2. A hadronic scenario 

We also considered a "hadronic scenario" for y-ray production 
from Cygnus X-3. In our model we used the same formalism 
adopted by Romero et al. (2003). In this case, the compact source 
is assumed to eject a flux of mildly relativistic hadrons (mostly 
protons) at the base of the jet. These protons are first acceler- 
ated near the compact object and then propagate along the jet 
interacting with the gaseous surroundings provided by the WR 
companion mass outflow. The resulting proton-proton (pp) col- 
lisions can copiously produce pions and y-rays resulting from 
neutral pion decays. 



14 Assuming that the photosphere radius is equal to the plasmoid ra- 
dius, i.e., 3 x 10 10 cm, the fraction of the plasmoid volume inside the 
y-ray photosphere is ~32%. 



The proton distribution in the jet is assumed to be isotropic 
in the jet comoving frame, with an energy spectrum described 
by a power law with a high energy cut-ofF: 

dN 

— — — Kp y exp(-y/y c ) [y > y m ,„] (2) 

dydv 

with a = 3, y„„„ = 1 and y c = 100. The spectral index of the 
distribution is the best-fitting value on the AG\LE-GRID spec- 
tral shape. We set the energy cut-off value at y c = 100 so that 
the total SED is consistent with the spectral constraints of the 
MAGIC upper limits. 

The ejected protons interact with the hadronic matter of the 
WR strong wind. The inelastic hadronic scatterings produce neu- 
tral pions that subsequently decay in y-rays. We adopted the 
same formula for the cross section cr pp (y) for inelastic pp in- 
teraction reported by Kelner et al. (2006). We assumed that the 
injected protons in the jet interact with the gas of the wind along 
a cylindrical column of matter 15 with a radius r — 3 x 10 10 cm 
and a height of H as 3 x 10 12 cm (this height provides the inter- 
esting part of the cylinder in which most of the interactions take 
place). In analogy with the leptonic models, we assumed for the 
jet a bulk Lorentz factor of Y — 1.5, and an inclination to the 
line of sight of i = 14°. In order to quantify the density of mat- 
ter in the WR wind, we assumed that the companion star has a 
mass loss rate of M ~ 10~ 5 M Q yr _1 and the speed of the wind 
is v wind ~ 1000 km s" 1 (Szostek & Zdziarski , 2008). By inte- 
grating the density of matter in this cylinder expressed in terms 
of number density of protons (g ~ 1 /R 2 , where R is the distance 
from the star), we find that the total number of protons from the 
wind in this column is N PtW i„d as 3.7 x 10 45 . 

We considered a variety of proton injection rates in the jet, 
Npjet- The result of our best hadronic model for Cygnus X-3 is 
reported in Figure 5. In this case, the integrated number of pro- 
tons injected in the jet is N p j et « 9. Ox 10 42 , and the average num- 
ber proton density in the column of interest is n p j et ~ 1.1 x 10 9 
protons cm 4 . 

The proton injection flux in the jet of our best model turns 
out to be <f> p jet ~ 2.4 x 10 19 protons ctrT 2 s , that corresponds 
to a proton injection rate of N p j e t ~ 6.7 x 10 40 protons s _1 . 

In Cygnus X-3 the corresponding jet kinetic luminosity for 
the hadrons (L^ p — N p T m p c 2 ) would be Z^,,, p w 1.5 x 10 38 
ergs -1 . This value is consistent with the average bolometric 
luminosity of the hypersoft state, L^ s as 1.2 x 10 38 erg s~ ! 
(Koljonen et al., 2010). Moreover, L,t,„ p is lower than the 
Eddington accretion limit for the system, that is L El td ~ 10 39 
erg s _1 , assuming that the compact object is a black hole with a 
mass of M x as 10M o . 

5. Discussion 

In the context of a leptonic scenario, we find that most of the 
y-ray emission above 100 MeV is due to IC scatterings on wind 
stellar photons (see Figures 3 and 4), according to the results of 
Dubus et al., 2010 and Zdziarski et al., 2012. We remark that the 
IC cooling times are very short (f /c ~1— 10 s). Thus, the observed 
time scale of y-ray emission (1-2 days) and the strong orbital 
modulation (~4.8 hours) detected in y-rays (Abdo et al., 2009) 
impose a continuous injection of accelerated particle in the jet. 
The y-ray emission, as remarked also by Zdziarski et al. (2012), 

15 A cylindrical configuration of the jet corresponds to set e = in the 
formalism adopted by Romero et al. (2003) to describe the jet radius 
dependence on the axis, r(z) = £,t . 
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Fig. 3: Multiwavelength SED of Cygnus X-3 during the main y-ray events (non-simultaneous data), leptonic model "A" (see main text). Blue 
circles: X-ray average "hypersoft" spectrum (Koljonen et al., 2010), RXTE-PCA and RXTE-HEXTE data (~3 to ~150 keV); red circles: AGILE- 
GRID energy spectrum (50 MeV to 3 GeV) of the main y-ray episodes (Figure A.l and A. 2); magenta arrows: MAGIC differential flux upper 
limits (95% C.L.), 199-3155 GeV, related to soft spectral state (Aleksic et al., 2010). Spectral components of the model: black body (BB) emission 
from disk (blue short-dashed line), Inverse Compton (IC) on soft photons of the accretion disk (green solid line), Inverse Compton (IC) on soft 
photons of the stellar wind (red solid line). Global SED model curve in black solid line. 
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Fig. 4: Multiwavelength SED of Cygnus X-3 during the main y-ray events (non-simultaneous data), leptonic model "B" (see main text). Spectral 
components of the model: BB emission from disk (blue short-dashed line), Inverse Compton (IC) on soft photons of the accretion disk (green 
solid line), Inverse Compton (IC) on soft photons of the stellar wind (red solid line). Global SED model curve in black solid line. For a detailed 
description of the datasets, see caption in Figure 3. 
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Fig. 5: Multiwavelength SED of Cygnus X-3 during y-ray events (non-simultaneous data), hadronic model (see main text). BB emission from disk 
(blue short-dashed line), y-ray emission from ^-"-decays (black solid line). For a detailed description of the datasets, see caption in Figure 3. 



cannot be related to a single injection of a relativistic plasmoid 
in the jet. Our leptonic picture could suggest a substructure. We 
find that the innermost part of the jet, where the density of X- 
ray disk photons is high, could give a strong contribution to the 
hard X-rays at ~ 100 keV. In model "A", the spectral component 
related to IC scatterings on soft photons from the disk (green 
curve) gives a substantial contribution to the overall model be- 
low ~100 MeV that is consistent with the hard tail of the X-ray 
dataset (see Figure 3). This model suggests a possible spectral 
link between the power-law tail of ultrasoft/hypersoft state and 
the y-ray emission detected by the AGILE-GRID. On the other 
hand, in the region of the jet outside the y-ray photosphere re- 
lated to pair production on disk photons (H > 10 10 cm), the "IC 
disk" component becomes very low at ~ 100 keV and the overall 
contribution to the hard X-rays is negligible (Figure 4). Thus, ac- 
cording to model "B", if the region of y-ray emission is far away 
from the compact object, the power-law tail in the hard X-ray 
band cannot be accounted for by IC processes in the jet. For this 
reason we assumed two "extreme" cases in the spatial configu- 
ration of the y-ray emitting blob, that is located at H « 3 x 10 10 
cm (model "A"), and at H » 3 x 10 12 cm (model "B"). Dubus 
et al. (2010) found that the y-ray emitting region is located at a 
distance H < lOd ~ 3 x 10 12 cm. By assuming that the compact 
object is a black hole with a mass of M x = 2OM , their simula- 
tion give the best fit of the y-ray modulation for H w 3 x 10 11 
cm. So that, the bulk of the y-ray emission is possibly produced 
at an intermediate configuration between models "A" and "B". 
The leptonic picture "B " is qualitatively consistent with the one 
proposed by Zdziarski et al. (2012) to explain the Fermi-LAT 
data published in Abdo et al. (2009), even if the peak energy of 
the IC bump is quite different: in our models the peak energy is 
between ~ 1 00 MeV and ~ 1 Ge V, in their models the peak energy 
is between ~3 and ~50 MeV. This difference could be due to 
the unequal spectral indices of the y-ray spectra detected by the 



AGILE-GRID and by Fermi-LAT (see Figure A.2). Moreover, 
Zdziarski et al. (2012) assumed that the electrons are injected in 
the jet with a power-law rate for j\ < y < j2- The electrons 
subsequently lose energy via Compton, synchrotron and adia- 
batic losses and form a distribution below j\ . Thus, they demon- 
strated that models with j\ < 10 3 give strong contribution to 
the hard X-rays, which appears in conflict with the observed or- 
bital modulation at ~ 100 kev during the y-ray emitting intervals, 
that is out of phase with the y-ray modulation at > 100 MeV 
(Zdziarski et al., in press). This phase misalignment would rule 
out any substantial contribution of the jet to the hard X-rays, in 
(apparent) contradiction with our findings in model "A". First, 
we remark that in our models the energy distributions of the ac- 
celerated particles are assumed to be "steady-state" spectra, sub- 
sequent to physical cooling processes, and not injected spectra. 
Then, our leptonic models are based on IC scatterings on both 
UV stellar photons and X-ray disk photons, whereas Zdziarski et 
al. (2012) consider the UV stellar photons only (neglecting the 
X-ray photons from the accretion disk). By observing Figure 3, 
we note that at ~ 100 keV the contribution of the "IC disk" com- 
ponent (green curve) is equivalent to the contribution of the "IC 
star" component (red curve). But, if the latter has a modulated 
emission (the jet-wind geometry is anisotropic during the or- 
bital phase), the former would not change with the orbital period 
(the jet-disk geometry should not significantly change with or- 
bital phase). Thus, the effective unabsorbed modulation at ~100 
keV - that is a superposition of a modulated and a unmodulated 
component - should have a lower amplitude with respect to the 
effective modulation at energies > 100 MeV - that is actually 
related to the modulated "IC star" component only. Moreover, 
anisotropic absorption effects in the wind could strongly affect 
the unabsorbed 100 keV emission and produce an observed mod- 
ulation in phase with the soft X-ray band. 
Finally, we remark that, for energies greater than ~10 GeV, y- 
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rays are above the threshold for pair production on stellar pho- 
tons. Nevertheless, Zdziarski et al., 2012 demonstrated that, for 
a similar choice of geometrical parameters, the value of the op- 
tical depth is moderate and peaks at ~0. 1—1 TeV. Thus, the yy 
absorption by UV stellar has been neglected in our (leptonic and 
hadronic) models. 

In the context of a hadronic scenario, we use a model similar 
to the one proposed by Romero et al. (2003). The only substan- 
tial difference consists in the jet geometry: we use a cylindri- 
cal model, whereas they use a conical configuration. We found 
that a simply hadronic model can account for the y-ray spec- 
trum detected by the AGILE-GRID, by assuming a reasonable 
proton injection rate in the jet. It is interesting to compare our 
best-fit hadronic injection rate for Cygnus X-3 with the value 
deduced for the microquasar SS433 that is known to produce 
a quasi-steady jet of hadronic nature (Migliari et al., 2002). 
SS443 is characterized by a hadronic jet mass ejection rates near 
M jel ~ 5 x 10~ 7 M Q yr l (Konigl, 1983; Fabrika & Borisov, 
1987; Reynoso et al., 2008), that corresponds to a proton injec- 
tion rate of ss 443 as 1.9 x 10 43 protons s _1 . Thus, we have 

that N Cyg - X ~ 3 * 3.5 x 10~ 3 N ss443 . SS433 ejects hadrons in a 

pjet pjet 

quasi-steady fashion, Cygnus X-3 ejects hadrons in highly vari- 
able regime with a lower injection rate. 

Our hypothesis for a hadronic interpretation of y-ray emission 
from Cygnus X-3 needs to be supported by information that 
at the moment is not yet available, such as hadronic emission 
lines in the flare spectra, and a precise characterization of the y- 
ray spectrum at energies below 100 MeV that should show the 
characteristic decrement of neutral pion emission. Furthermore, 
hadronic mechanisms, besides emitting strong y-ray radiation 
via 7r°-decay, would produce intense flux of high-energy neu- 
trinos, emerging from the decay of secondary charged mesons 
produced in pp collisions. So, a firm simultaneous detection of 
strong neutrino flux and y-ray activity from Cygnus X-3, would 
represent the signature of a dominant hadronic mechanism in the 
relativistic jet. 

In this scenario, considering the temporal coincidence of y- 
ray/radio flares, we are implicitly assuming that the hadronic 
component of the jet provides the main contribution to the y- 
ray emission, and the leptonic component produces - via syn- 
chrotron emission process - the strong radio flares far away 
from the compact object. In addition to y-rays from ^"-decays, 
hadronic pp interactions are expected to produce a population of 
secondary electrons (and positrons) from the decay of charged 
pions (tt*). These secondary leptons can contribute to the emis- 
sion in the radio band via synchrotron processes and in y-rays 
(marginally with respect to the contribution by 7r° -decay) via IC 
and Bremsstrahlung processes. 

6. Conclusions 

Several events of y-ray activity were detected by the AGILE- 
GRID from Cygnus X-3 while the system was in a special 
radio/X-ray spectral state: intense y-ray activity was detected al- 
ways during prominent minima of the hard X-ray light curve 
(corresponding to strong soft X-ray emission), a few days be- 
fore intense radio outbursts (major radio flares). This tempo- 
ral repetitive coincidence turns out to be the spectral signature 
of the y-ray activity from this puzzling microquasar, and might 
open new areas to study the interplay between the accretion disk, 
the corona and the formation of relativistic jets. The simultane- 
ous strong soft X-ray emission from the disk and y-ray emission 
from the jet preceding intense radio outbursts might suggest a 



scenario in which the hot thermal corona "dissolves", and the 
accretion power from the disk directly charges the jet, that then 
emits y-rays and, subsequently, radio outbursts (via synchrotron 
processes) far from the compact object. 

The y-ray detections of Cygnus X-3 imply new constraints 
to the emission models for such powerful X-ray binary, indi- 
cating that hybrid-Comptonization mechanisms (Coppi, 1999) 
alone cannot account for the y-ray fluxes detected by AGILE 
and Fermi above 100 MeV, unless we assume unrealistic physi- 
cal parameters (Cerutti et al., 201 1). This implies that the corona 
cannot be the site of y-ray emission. We found that the innermost 
part of the jet (distances < 10 10 cm from the compact object) 
could give a strong contribution to the hard X-rays at ~ 100 keV 
during the y-ray emitting interval, while the farthest part (dis- 
tances > 10 1() cm from the compact object) produces the bulk of 
the y-ray emission above 100 MeV. 

We found that the y-ray spectrum of Cygnus X-3 detected 
by the AGILE-GRID is significantly harder than the time aver- 
aged spectrum obtained by Fermi-LAT for the "y-ray active pe- 
riods" of the microquasar, lasting ~4 months (see Figure A. 2). 
Although both the AGILE main y-ray events and the Fermi y- 
ray active periods are both likely related to the presence of an 
active jet, the spectral difference may imply a fast hardening of 
the spectrum during the peak y-ray events, lasting ~l-2 days. 

We have demonstrated that both a leptonic model based on 
inverse Compton emission from a relativistic plasmoid injected 
in the jet and a hadronic model based on 7r°-decays, might ac- 
count for the y-ray emission observed by the AGILE-GRID. 
Both these models require the introduction of a new component 
("IC bump" or "7r°-bump") in the SED of the system. 
A leptonic scenario seems to be more likely than a hadronic one: 
the y-ray modulation, the spectral link between hard X-ray and 
y-ray spectra, the temporal link between y-ray events and ra- 
dio flares could be interpreted in a natural way by assuming that 
the electrons are the main emitters. According to our results, the 
y-ray emission occurs within ~10 12 cm from the compact ob- 
ject. Interpreting as the propagation time of the relativistic jet 
(v = V5c/3) the ~4-day delay between the onset of y-ray and 
radio flaring emission, the radio burst would occur at a distance 
of ~8 x 10 15 cm. 

Our hadronic model, with the assumption of a standard WR 
wind, would require a jet kinetic power of p m 1.5 x 10 38 
erg s _1 to explain the y-ray emission detected by AGILE. This 
value is of the same order of magnitude as the bolometric lumi- 
nosity of the disk/corona during the hypersoft spectral state, and 
it is lower than the Eddington accretion limit for a black hole 
with a mass of M x * 10M o (L Edd * 10 39 erg s -1 ). So that, a 
hadronic picture is physically reasonable and not energetically 
disfavored with respect to a leptonic one. At present, there is no 
strong evidence for excluding one of these hypotheses, and it 
remains an open question whether the dominant process for y- 
ray emission in microquasars is hadronic or leptonic (Mirabel, 
2012). 

The firm discovery of y-ray emission from this microquasar 
represents the experimental proof that these astrophysical ob- 
jects are capable of accelerating particles up to relativistic ener- 
gies, through a mechanism - related to the disk-corona dynamics 
- that leads to jet formation. 
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Appendix A: The AGILE-GR/D dataset 

We performed an analysis on the whole AGILE-GRID data in 
the period November 2007 - July 2009, using a detection al- 
gorithm that was developed by the AGILE Team to automati- 
cally search for transient y-ray emission. The algorithm initially 
analyzed 140 maps, each related to a 2-day integration (non- 
overlapping consecutive time intervals). The time bins of peak y- 
ray emission, with significance of detection greater than 3 sigma, 
have been found. The analysis was subsequently manually re- 
fined in order to optimize the determination of the time interval 
of the y-ray emission. The whole analysis was carried out with 
the Build 19 version of the AGILE team software, using the 
FM3 . 119_2 calibrated filter applied on the consolidated dataset 
with off-axis angles less than 40°. We used a Multi-Source 
maximum-Likelihood Analysis (MSLA) to take into account the 
emission of the nearby y-ray pulsars 16 : 1AGL J202 1+3652 (PSR 
J2021+3651), 1AGL J2022+4032 (PSR J2021+4026), 1AGL 
J2032+4102 (PSR J2032+4127). In particular the MSLA is fun- 
damental in order to avoid contamination by the pulsar PSR 
J2032+4127, located at a distance of ~ 0.5°. In this paper 
we did not consider the off-pulse data for the nearby pulsar. 
Nevertheless, the MSLA accounted for the steady y-ray emis- 
sion from this source in the calculation of the significance and 
the flux of each y-ray detection of Cygnus X-3. Moreover, we 
can exclude any substantial spectral contamination from the pul- 
sar because the steady y-ray emission from the pulsar 17 is much 
smaller than the mean flux of the active y-ray emission from 
Cygnus X-3. 

The main events of y-ray activity, detected with a significance 
above 3<x ( y/TS > 3), are shown in Table A.l. We found seven 
events 18 , including those presented in Tavani et al. (2009a) and 
in Bulgarelli et al. (2012a). 

By integrating all the main events with FM3 . 1 19_2 filter, we 
found a y-ray source detected at 6.7cr ( y/TS = 6.7) at the aver- 
age Galactic coordinate (/, b) = (19.1°, 0.9°) ± 0.4° (stat) ± 0.1° 
(syst), with a photon flux of (158 + 29) X 10" 8 photons cnT 2 s _I 
above 100 MeV 19 . The average differential spectrum between 
100 MeV and 3 GeV is well fitted by a power law with a photon 
index a = 2.0 + 0.2 (Figure A.l). This value is consistent with 
the Cygnus X-3 photon index found by Tavani et al. (2009a) and 
Bulgarelli et al. (2012a). In Figure A. 2 we present a compari- 
son between the vF v spectra of Cygnus X-3 as detected by the 
AGILE-GRID and by Fermi-LAT (Abdo et al., 2009) during the 
y-ray activity. It should be noted that the AGILE-GRID spec- 
trum is related to the peak y-ray activity only (the seven main 
events, lasting 1-2 days, in Table A.l), instead the Fermi-LAT 



16 The main characteristics of the persistent y-ray sources we used in 
the MSLA are reported in the Table 1 of Chen, Piano et al. (201 1) 

17 The steady y-ray emission from the pulsar PSR J2032+4127, as 
detected by the AGILE-GRID, is F* ss = [37 ± A(stat) ± I0%(syst)] x 
10~ 8 photons cirr 2 s -1 for photon energy above 100 MeV, see Chen, 
Piano et al. (201 1) for details. 

18 All the y-ray events were detected by using the same filter, 
FM3.119_2. As discussed in the Supplementary Information of Tavani 
et al. (2009a), the event of 2-3 November 2008 appears to be relatively 
"soft" in y-rays compared to the other episodes. In fact, by analyzing the 
event with the FT3ab_2 filter (more efficient in detecting such kind of 
emission), we found a more significant detection: y/TS = 3.9, photon 
flux above 100 MeV: (214 ± 73) x 10~ 8 photons cnr 2 s _1 . 

19 Here we present an updated result for the analysis of the 7-events 
integration with respect to the one reported in Piano et al. (2011). In 
this paper the analysis has been carried out with a more recent version 
of the AGILE software tool (AG_multi4). 



spectrum is an average spectrum detected during the two active 
windows (2 months each) of y-ray emission from Cygnus X-3 
(MJD: 54750-54820 and MJD: 54990-55045). 




Energy [MeV] 



Fig. A.l: Photon spectrum between 50 MeV and 3 GeV of Cygnus X-3 
as detected by the AGILE-GRID by integrating all the main y-ray 
episodes in Table A.l. Power law fit on y-ray data between 100 MeV 
and 3 GeV with photon index a = 2.0 ± 0.2. 



-9,0 n 




log [E (MeV)] 

Fig. A.2: vF y spectra of Cygnus X-3 during the y-ray activity. Red cir- 
cles: AGILE-GRID energy spectrum (50 MeV to 3 GeV) of the main 
episodes (Figure A.l). Blue error contours and dashed blue line: av- 
erage power-law fit with a = 2.70 ± 0.25 of the spectrum detected by 
Fermi-LAT during long integration (~4 months) in the active periods 
(Abdo et al., 2009). 



Finally, we evaluated the post-trial significance for repeated 
flare occurrences by using the same formalism of Bulgarelli et 
al. (2012b). The probability of having k or more detections - 
consistent with the position of Cygnus X-3 - with ^JTS > yfh 
in TV trials, is: 

where p is the p-value corresponding to h. For y/FS > 3 we 
have a p-value of p = 2.0 x 10~ 3 . So, for N - 140 (our tri- 
als, the number of 2-day integration maps analyzed by the initial 
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Table A.l: Main events of y-ray emission detected by the AGILE-GRID in the period November 2007 - July 2009. All detections have a signif- 
icance above 3<r ( s/TS 3> 3). Column one: period of detection in UTC; Column two: period of detection in MJD; Column three: significance of 
detection; Column four: photon flux (above 100 MeV). 



Period 


MJD 


Vts 


Flux [10 8 photons cm 2 s '] 


2008-02-11 (18:07:28) - 2008-02-12 (11:07:44) 


54507.76 - 54508.46 


3.7 


264 ± 104 


2008-04-16 (13:59:12) - 2008-04-17 (13:48:00) 


54572.58 - 54573.58 


4.5 


265 ± 80 


2008-11-02 (13:01:05) - 2008-11-03 (19:01:05) 


54772.54 - 54773.79 


3.1 


135 ±56 


2008-12-11 (19:50:40) - 2008-12-12 (23:02:40) 


54811.83 - 54812.96 


4.0 


190 ± 65 


2009-06-20 (21:04:48) - 2009-06-21 (20:53:04) 


55002.88 - 55003.87 


3.8 


193 ± 67 


2009-07-13 (01:11:60) - 2009-07-14 (00:59:44) 


55025.05 - 55026.04 


3.2 


216 ± 89 


2009-07-21 (21:07:12) - 2009-07-23 (21:07:12) 


55033.88 - 55035.88 


3.6 


158 ±59 



algorithm) and k — 7 (our detections of Cygnus X-3), we found 
P(140,7) = 1.8xl0~ 8 , that corresponds to 5.5 Gaussian standard 
deviations. 

A Multi-Source maximum-Likelihood Analysis (MSLA), 
applied in the deep integration of the AGILE-GRID data (be- 
tween November 2007 and July 2009), found a weak persistent 
emission (significance: y/TS = 5.17) from a position consistent 
with Cygnus X-3 20 with a photon flux of (14 + 3) x 10~ 8 photons 

-2 -1 

cm s . 
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